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Abstract. We discuss recent progress in the study of in-medium heavy quarkonium using 
first-principles lattice QCD calculations and effective field theory. In particular our focus lies 
on real-time information carried by QCD spectral functions and we report on a novel strategy 
for their extraction from Euclidean-time simulations. Combining an effective description for 
heavy quarks with finite temperature lattice QCD we compute in-medium spectra providing 
insight in heavy quarkonium melting in a static thermal medium. On the other hand we use 
spectral information to extract the complex in-medium potential which provides the basis for a 
dynamical description of quarkonium real-time evolution via the Schrodinger equation. 


1. Introduction 

Heavy quarkonium, the bound state of a b or c-quark with its anti-partner, is a valuable probe 
of the strong interactions both in vacuum, as well as at high energy densities, such as in the 
droplet of quark-gluon-plasma created in relativistic heavy-ion collisions at RHIC and LHC [T]. 

Recent experimental studies of vector channel Bottomonium (T(reS')) and Charmonium 
(J/'L,'l'(nS')) by the CMS and ALICE collaborations have revealed two distinct faces of 
quarkonium, represented by the different flavors each. For bb, a comparison of di-muon spectra 
inp+p and Pb + Pb collisions provided evidence for a pattern of sequential suppression, i.e. less 
deeply bound excited states are dissolved more easily than the ground state [2j. These findings 
are in line with the concept of Bottomonium forming early on in a collision and traversing the 
plasma as well defined probe. On the other hand ALICE observed [3] at mid-rapidity that the 
suppression of Charmonium reduces with increasing centrality. This supports the picture of J/T 
formation from fully thermalized cc at the phase boundary, as advocated for in the statistical 
model of hadronization [3]. Understanding the intricate in-medium physics leading to each of 
the two scenarios from first principles has to be the central goal for theory. 

Heavy quarkonium is characterized by the well pronounced separation of scales between the 
constituent rest mass {rric = 1.275(25)GeV, = 4.66(3)GeV) [S] and all other scales present in 

current experiments. This distinguishes them as well defined observable, measurable with high 
precision, and opens up the possibility for theory to employ effective descriptions that allow a 
more intuitive understanding of the bound state physics in the presence of a thermal medium. 

Already at zero temperature it is the intrinsic scale of quantum fluctuations of QGD, often 
referred to as Aqcd — 200MeV that lies well below mg, which in turn tells us that quantum 
fluctuations cannot spontaneously create a QQ pair. Such a system, where creation and 


annihilation do not play a central role can be described by a non-relativistic effective field- 
theory (EFT) called NRQCD, in fact one can show that it can even be described to a good 
degree by a non-relativistic potential entering a Schrodinger equation . At finite temperature 
the additional scale T enters and one might ask how its presence changes the properties of the 
bound state, which eventually will melt if thermal fluctuations overcome its typical binding 
energy. NRQCD when used at T=0 only requires a separation of tuq from the next lower scale 
and thus in principle is valid also in-medium. We report on its application to lattice QCD studies 
of heavy quarkonium at T > 0 in more detail in secj^ 

On the other hand, the question whether an in-medium potential V{r) can be defined and 
what its values are, has only recently been answered in a satisfactory way. For more than two 
decades, in the absence of the means to derive a potential from QCD, theorists have resorted 
to model potentials, among which two candidates have gained popularity. The purely real 
color singlet free energies in Coulomb gauge = (Tr D(x)Dl(0) )cg from Polyakov-loop 


D(x) = exp ig Jq drA'^^Tjx) correlators [7j and the internal energies 
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from it [8], as well as linear combinations [9] have been employed. 

It is important to note that while these quantities are readily accessible in lattice QCD |10j 
they have not been linked to a potential in QCD. In fact progress in the formulation of EFT’s 
mm made it possible to actually derive an expression for the potential between two static 
quarks in a thermal medium starting from the QCD Lagrangian, revealing that neither nor 
alone can play the role of V (r). 

The reason is that the evaluation of the potential both in resummed hard-thermal loop 
perturbation theory [laii! and on the lattice |15[ fT6] showed that V (r) is actually a complex 
quantity. While at T > Tc, Re[R] is closely related to the phenomenon of Debye screening, Im[I/] 
is intimately linked to the scattering of the medium partons with the color string spanning 
in between the QQ m m- Using a purely real model potential neglects these effects and 
necessarily misses an important aspect of the in-medium physics of heavy quarkonium. Recent 
lattice QCD results on the extraction of the complex in-medium potential are presented in sec|^ 
While the non-relativistic character of the heavy quarks allows a simplification of their 
description, the confinement aspects of the QCD vacuum as, well as the physics of the strongly 
coupled QGP at T ~ Tc require fully relativistic, as well as non-perturbative methods. One 
truly first principles approach is lattice QCD, where spacetime is discretized on a hyper-cubic 
grid with finite lattice spacing a and box size N^a. Gauge fields denoted in continuum by 
are represented by SU(3) valued link variables U^{x) = exp[ig'aA“T“], where T“ represent the 
generators of the gauge group. These object reside on the links connecting the nodes of the 
spacetime grid. The light fermions of the QGD medium on the other hand populate the nodes. 

The strength of the lattice approach lies in the fact that the Feynman path integral of this 
theory can be evaluated non-perturbatively via stochastic Mont-Garlo methods, after rotating 
the Minkowski-time axis onto imaginary time. Since the simulation is not carried out in physical 
time a form of analytic continuation is required if dynamical information is to be extracted. This 
represents an inherently ill-defined problem. One strategy, Bayesian inference, proven successful 
in many areas of physics to overcome such a challenge is discussed in secj^ 

At T = 0 the combination of lattice QGD and NRQGD, have been highly successful in 
reproducing the quarkonium spectrum (see e.g. [13 [E]). There is no modeling input necessary in 
these computations, thus they truly link the QGD Lagrangian to physically observable quantities, 
such as particle masses or matrix elements. Note that since all simulations are carried out 
using dimensionless numbers, at some point experimental input is required to set the overall 
scale. Once the scale is fixed however ab-initio predictions are possible, one example being the 
existence of the %(2S) state with mass ^ = 9.988(3)GeV predicted in 2011 [IB], which 

was confirmed by the Belle collaboration [19] at ~ 9.999(3.5)(2.8) in the following year. 






At finite temperature lattice QCD can also help us to elucidate the static properties of 
the QCD medium, now that consistent continuum extrapolated results for thermodynamic 
quantities, such as the energy density and pressure have been obtained [HI Ea [22]. While 
the transition between the confined and deconfined phase is a mere crossover at physical 
quark masses m, a pseudo-critical temperature of Tq = 155(9)MeV has been established that 
characterizes the change from the hadronic to QGP regime. Evaluation of the trace anomaly 
[saiiaisii, a non-perturbative measure of the interaction strength revealed a region of strong 
coupling exists in the regime just above Tq- This reinforces the need for non-perturbative 
methods when approaching the physics of heavy quarkonium in heavy-ion collisions. 

2. Spectral functions from lattice QCD 

Physical information relevant for an understanding of in-medium heavy quarkonium can be found 
in spectral functions p{uj). They inform us about the presence and strength of bound states, 
which appear as peaked features of skewed Breit-Wigner type. The position encodes the mass of 
the state, the width its transition probability to other bound states or into the continuum, which 
is manifested as a featureless structure at large w. What makes spectral functions interesting 
from a technical point of view is that they constitute a bridge between correlation functions in 
real-time and those analytically continued to Euclidean time, which will play a crucial role in 
the extraction of the heavy quark potential from lattice QCD simulations in secj^ 

In secj^we will calculate heavy quarkonium correlation functions D{t) on the lattice. These 
are related to spectra via a Kernel which is determined by QCD and exactly known 

p ^UJ 

D{t,x)= duj K{t,uj) p{u},x), D{Ti) = AuiiKii pi. (1) 

1=1 

In both studies presented in secj^ and sec|^ we will encounter the T independent Kernel 
K(t,uj) = , with which the above expressions amount to Laplace transforms. 

Writing the convolution on the right of Eq. 0 with Nr steps along imaginary time and Ni^ 
frequency points, reveals the challenge we face in inverting this relation. I.e. a function p with 
an intricate set of features has to be computed from a finite set of noisy datapoints obtained from 
a Monte Carlo simulation. This is an ill-defined problem, since for any discretization r* = /3/Al¬ 
and in the presence of inevitable noise on D, an infinite number of choices for the p/s exist that 
will reproduce the data, even in the presence of the unitarity constraint p; > 0. 

One possibility to give meaning to the problem is through so called Bayesian inference, which 
is based on Bayes theorem P[p\D,I] oc P[D|p,/]P[p|I]. It tells us that the probability for a 
test function p to be the correct spectral function, given the simulated lattice data D and any 
further prior information I on the spectrum, consists of essentially two terms. The likelihood 
probability P[D\p,I] = e~^ denotes the probability of the lattice data if p were the correct 
spectrum and L amounts to nothing but the usual quadratic distance functional. The second 
term -P[p|/] = denotes the prior probability, which encodes prior information on the spectrum 
and acts as regulator of the degenerate fit- On the one hand the functional form of S encodes 
part of I, on the other hand ^[p, m] depends on a function m{uj), the default model, which 
defines the extremum of S and thus represents the correct spectrum in the absence of data. We 
have to carry out a numerical optimization to find the most probable p given D and /, with the 
competing L and S leading to a unique answer. 

One popular implementation of the Bayesian strategy is the Maximum Entropy Method 
(MEM) [25l |26l |27|. In it the Shannon-Jaynes entropy is used as regulator term S = 5sj. 
One difficulty lies in the fact that S'sj contains flat directions, which can adversely affect the 
convergence properties of the method and in practice necessitates the use of a priori restrictions of 
the functional space from which the spectra can be chosen. To avoid this resolution impediment 
we instead use a recently developed Bayesian approach [28j that relies on a prior S derived from 


the following three physical consideration: positivity; smoothness of the spectrum, where data 
has not introduced sharp features and the independence of the result from the units chosen. As 

) > which 


discussed in [28] these lead us to a prior functional S = a ^ ^ 


mi 


is devoid of the flat directions of 5sj. Note that the weighting factor a here in contrast to the 
MEM can actually be integrated out explicitly in a Bayesian fashion |28] . With S set, we only 
need to specify m{oj), for which we will use a constant, the most neutral choice. 


3. Bottomonium spectra from lattice NRQCD 

Let me paint in broad strokes the ingredients to investigating heavy-quarkonium in the non- 
relativistic field theory NRQCD on the lattice at finite temperature using Bottomonium as 
example. Readers interested in more technical details can consult e.g. m- Due to the 
separation of scales between heavy-quark rest mass and the two other characteristic scales Aqcd 
and T, creation and annihilation of heavy quark pairs is highly suppressed. This allows us to 
go over from a description of the quark fields in terms of four-component Dirac spinors Q{x) 
to a set of two two-component Pauli spinors representing the quark and anti¬ 

quark contribution respectively. Formally this separation of the upper and lower components 
of the four-spinor is obtained from a systematic expansion of the QCD Lagrangian in terms of 
increasing powers of the inverse rest mass (or more precisely the heavy quark velocity v) via the 
Foldy-Tani-Wouthuysen transformation. This is the same approach used to derive the Hydrogen 
energy levels from the Dirac equation. 

If we place the non-relativistic theory on a Euclidean spacetime lattice, the Pauli spinors 
obey the corresponding non-relativistic Hamiltonian, 

H = Ho + 6H = - —cr • B + -^(A^ • E - E • A±) + ... (2) 

2Mb 2Mb 8M2 ^ ^ ^ ^ 

Here and and denote the lattice derivative and Laplace operator respectively. If all 
factors of the lattice spacing a were made explicit it turns out that the expansion parameter 
amounts to aMb on the lattice. For the UV behavior to be well defined its value must be larger 
than 1.5 [29], which reflects the fact that NRQCD works better on coarser lattices. 

The fact that quark and antiquark modes are separately propagating in the medium 
background leads to an initial-value problem for the heavy quark propagator G(x, rQ = 
(V’(x, rj)'!/>(0, 0)) in the background of medium fields. Written in terms of the link variables 
provided as external fields by conventional lattice simulations we have 

G(x,ri) = H](x,r)^l-^^ x (1 - 5iL) G(x, ri_i) 

The parameter n reflects a freedom in the choice of the time discretization and is set to 
n = 2 here. From G(x, Tj) we can construct heavy quarkonium propagators D(x, r) = 
(0(x, r)G(x, r)Ol(xo, ro)G^^(x, r))med selecting the individual channels by appropriate 

O O 

vertex operators 0('*5i;x, r) = ai for the S-wave 0('^Pi;x, r) = AjCJj — AjCJj for the P-wave. 

The medium expectation value is obtained by measuring and averaging these correlators 
on realistic gauge field configurations containing light u,d, and s (HISQ) quarks, generated 
by the HotQCD collaboration [2T|. These 48^ x 12 lattices with Tc = 159MeR feature 
an almost physical pion mass of mn = 161MeV and allow a fine temperature scan between 
T = 140 — 249MeV, since the lattice spacing is used to regulate the extend in Euclidean time 
/3 = 6.664 — 7.280. The expansion parameter lies between aMb = 2.76 — 1.56 > 1.5 allowing for 
the application of NRQCD. For a series of studies using a complementary fixed scale approach 
see [30]. Experimental input, e.g. the T(IS') mass is required to set the relative and absolute 
frequency scale, for which we use low temperature ensembles with 48^ x 48, 64 . 
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Figure 1. S-wave (left) and P-wave (right) Bottomonium spectra at T ~ 0. 
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Figure 2. Spectra of the full (solid) and truncated (dashed) T = 0 correlator sets. 


At T PS 0 we have extracted [32] (left) the S- and (right) P-wave spectra shown in Figj^ 
The S-wave ground state peak is very well resolved, around an order of magnitude more clearly 
than what is possible with the MEM. The second bump contains mostly contributions from 
T(25), the higher lying wiggles are numerical ringing. After fixing the scale by the S-wave, 
we can estimate the systematic error in the computation from our postdiction for the Xb(lS') 
mass. We obtain a value of = 9.917(3)GeV slightly larger than the PDG value of 

m™^p) = 9.89278(26)(31)GeV. The P-wave ground state appears more washed out than the 
S-wave, not only because of a lower signal to noise ratio in the underlying correlators. Even in 
continuum the ratio of P-wave peak height to background is smaller in comparison m- 

To investigate in-medium modification we have to ascertain how well our spectral 
reconstruction fares at T > 0. The difference to vacuum lattices Nj- = 48 is the physical 
extend in Euclidean time, reflected in a smaller number of Nr = 12 correlator points at 
r > 0. Hence we take the T = 0 correlator data set and truncate it after twelve points 
{K{uj, t) is T independent). Supplying it to the reconstruction yields the spectra shown as gray 
dashed curves in Figj^ On the coarsest lattices (/? = 6.664) where NRQCD works best, the 
ground state peak reconstruction remains unaffected (Am < 2MeV,Ar < 2MeV), while on the 
finest lattices (/3 = 7.280) corresponding to the highest temperature in the following, we find 
(Am < 40MeV,Ar < 21MeV). As expected for the P-wave the accuracy limits are even larger 
[32] . Note that the changes are on the level of individual MeV, which could not have been 
resolved in the same manner with the MEM. If at finite temperature we were to find changes 
larger than the stringent limits set here, we can attribute them to in-medium modification. 

In Figj^ we show the T > 0 resnlts for (left) the S-wave and (right) the P-wave spectra at 
the fourteen different temperatures available. In the S-wave we find at all T a well defined 
ground state peak, much larger than any further numerical ringing at higher u. If zoomed in, 
the peak is seen to slightly shift and broaden with T, however the changes are smaller than the 
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Figure 3. S-wave (left) and P-wave (right) Bottomonium spectra at T > 0. 


accuracy limits established before. Hence we refrain from attributing them to medium effects. 
Nevertheless by a simple inspection by eye it is evident that T(IS') survives up to T = 249MeV. 
This conclusion agrees with the results from recent a MEM based study m- 

The P-wave case is more subtle, not only because of the larger errors due to the lower signal 
to noise ratio in the correlators. While we find a well defined lowest lying peak at all T its height 
is comparable to that of higher lying numerical ringing at the highest T. A naive inspection 
by eye fails to give conclusive insight. Therefore we propose to put the decision of survival or 
melting on more systematic footing via a comparison with non-interacting spectral functions. 

The QQ spectrum in the absence of 
interactions is analytically known, not to show 
any bound state peaks. Such spectra can also 
be extracted numerically from lattices with all 
links set to unity and subsequent Bayesian 
reconstruction. Just as with Gibbs ringing 
in the Fourier transform we expect wiggly 
structures to appear, even if absent in the 
actual free spectrum. This can be seen in the 
gray curve in Figj^ where the free spectrum 
has been shifted to match the uj range of 
the full spectrum (red). We see that while 
numerical ringing dominates the second bump, the first peak is a factor of three larger and thus 
constitutes a signal for bound state survival. Hence we conclude that Xbi^P) survives up to 
T = 249MeV (l.GlTc). This conclusion differs from a recent MEM based study, which found 
melting already at T = 1.27Tc [30]. One explanation might lie in the restricted search space of 
the MEM, which shrinks with the number of available datapoints and thus entangles physical 
effects and numerical limitations. 
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Figure 4. Comparison of free (gray), 
interacting (red) and the spectral function from 
the truncated T = 0 dataset (blue dashed) 


4. The complex in-medium heavy quark potential from lattice QCD 

We now wish to address the question of how to determine the potential V (r) acting between two 
static quarks at finite temperature non-perturbatively. One starts from an EFT description in 
terms of QQ wave-functions called pNRQCD mm- The values of V (r) are determined by the 
procedure of matching, in which a suitable QCD correlation function is identified, which carries 
the same physics content as a correlation functions in pNRQCD and both quantities are set equal 
at a certain scale. At mg —>■ oo the correlator of two color-singlet wavefunctions is identified 
with the real-time Wilson loop Wn{t, r). The potential arises from its late Minkowski-time limit 

Wn{t,r) = {TT(exp[ig V{r) = ^lim 


(3) 

















Figure 5. (left) Real- and (left) imaginary part of the complex in-medium potential extracted 
from quenched QCD simulations. 


The limit reflects the fact that replacing an interaction mediated by gluons with an instantaneous 
potential is only valid at timescales much larger than those of the medium gluons. 

The challenge to a non-perturbatively evaluation lies in the real-time character of the 
definition, not directly amenable to a lattice QCD simulation. Here the concept of spectral 
function comes to our rescue |33l l3^ since the Wilson loop can be rewritten as a Fourier 
transform of a positive dehnite /9n(a;,r). In fact the Euclidean Wilson loop is related to the 
same pn via the Laplace transform, which can be inverted using the Bayesian method of 
Inserting the spectral decomposition in the potential dehnition ([^ leads to V{r) = 


sec 


limt^oo / doj p{u],r)/ f duj e~^‘^^p(u!,r). The late time limit tells us that only the lowest 

lying spectral structure contributes to V{r), which has been shown to take on the form of a 
skewed Lorentzian m- Its position and width correspond to the real- and imaginary part of 
the potential respectively. Hence we fit the lowest lying peak and read off Re[F] and Im[F]. 

In Figj^we show (colored points) the results |15[I16| for (left) the real- and (right) imaginary 
part obtained from anisotropic 32^ x Nj- lattices with /3 = 7 and bare anisotropy ^ = 3.5 in 
the quenched approximation (Tc = 270MeV). Changing W = 24 — 96 we scan a range of 
T = 210 — 839MeV. Note that instead of Wilson loops we here use Wilson line correlators in 
Coulomb gauge since they do not posses cusp divergences. 

Re[F] shows a smooth transition from the well 
known Cornell type behavior at low T to a Debye 
screened form at T > Tc, as expected from 
perturbation theory. Interestingly the values of (r) 

(gray points) lie very close to the values of the real part. 

Hence this first principles calculation supports 
as an approximation of Re[F]. Im[F] is non-zero on 
the lattice at T > Tc and hence cannot be neglected 
in an in-medium potential description. Its shooting off 
at larger r is an artifact of the reconstruction due to a 
diminishing signal to noise ratio in the correlators. 

We present in FigJ^ the values of Re[I7] in a true 
QGP with u,d, and s quarks. It is extracted on isotropic 
48^ X 12 lattices by the HotQCD collaboration [T^ [T6] . 
spanning a temperature range of T = 148 — 286MeV corresponding to a lattice parameter range 
of /? = 6.8 — 7.48. Based on the asqtad action with ~ 300MeV the transition temperature 
here is T = 172.5MeV. The behavior we find is qualitatively similar to the quenched case with 
a smooth transition from Cornell to Debye screened behavior, and again F^^'^ (r) appears to be 
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Figure 6. Re[I7] in full QCD with 
dynamical u,d and s quarks. 



































a good approximation of Re[y]. We refrain from showing the tentative values of Im[l/], since 
spectral widths extracted from W = 12 points are generally considered not reliable. 


5. Conclusion 

Heavy quarkonium presents a formidable testing ground for our capability to connect a hrst 
principles understanding of the strong interactions to phenomenologically relevant observables. 
The combination of non-relativistic EFT’s with non-perturbative simulations of the QCD 
medium, already highly successful at T ~ 0, are nowadays providing valuable insight also at 
T ~ Tq, reached in relativistic heavy ion-collisions. We reported on recent progress made in the 
determination of dynamical information from such simulations through a novel Bayesian spectral 
reconstruction approach and two studies based on its application were discussed. The hrst 
work extracted the in-medium spectra of bottomonium directly from lattice NRQCD and found 
survival of both S-wave and P-wave ground states up to T = 249MeV. The second work used 
spectral functions to extract the values of the complex in-medium heavy-quark potential. Besides 
a smooth transition from a conhning to Debye screened behavior it conhrmed the existence of 
an imaginary part in the potential in the quark-gluon-plasma. 

The author thanks Y. Burnier, O. Kaczmarek, S. Kim and P. Petreczky for the fruitful 
collaboration that underlies the lattice QCD spectra results presented in this proceeding. 
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